ApJ, IN PRESS; ASTRO-PH/0402247 

Preprint typeset using I4TgX style emulateapj v. 2/10/04 



O 
O 

(N 

X> 
PL. 



> 

^1" 

(N 
(N 
O 

o 

Or 

o ■ 



13 



CHANDRA OBSERVATIONS OF THE "DARK" MOON AND GEOCORONAL SOLAR- WIND CHARGE TRANSFER 

B. J. Wargelin 1 , M. Markevitch 2 , M. Juda, V. Kharchenko, R. Edgar, and A. Dalgarno 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
ApJ, in press; astro-ph/0402247 

ABSTRACT 

We have analyzed data from two sets of calibration observations of the Moon made by the Chandra X-Ray 
Observatory. In addition to obtaining a spectrum of the bright side that shows several distinct fluorescence 
lines, we also clearly detect time-variable soft X-ray emission, primarily O VII Ka and O VIII Lya, when 
viewing the optically dark side. The apparent dark-side brightness varied by at least an order of magnitude, up 
to ~ 2 x 10~ 6 phot s" 1 arcmin" 2 cm" 2 between 500 and 900 eV, which is comparable to the typical 3/4-keV-band 
background emission measured in the ROSAT All-Sky Survey. The spectrum is also very similar to background 
spectra recorded by Chandra in low or moderate-brightness regions of the sky. Over a decade ago, ROSAT also 
detected soft X-rays from the dark side of the Moon, which were tentatively ascribed to continuum emission 
from energetic solar wind electrons impacting the lunar surface. The Chandra observations, however, with their 
better spectral resolution, combined with contemporaneous measurements of solar-wind parameters, strongly 
favor charge transfer between highly charged solar-wind ions and neutral hydrogen in the Earth's geocorona 
as the mechanism for this emission. We present a theoretical model of geocoronal emission and show that 
predicted spectra and intensities match the Chandra observations very well. We also model the closely related 
process of heliospheric charge transfer and estimate that the total charge transfer flux observed from Earth 
amounts to a significant fraction of the soft X-ray background, particularly in the ROSAT 3/4-keV band. 
Subject headings: atomic processes — Moon — solar wind — X-rays: diffuse background — X-rays: general 



1. INTRODUCTION 

As reported bv lSchmitt et al.l 1199 lb . an image of the Moon 
in soft X-rays (0.1-2 keV) was obtained by the Rbntgen 
Satellite (ROSAT) using its Position-Sensitive Proportional 
Counter (PSPC) on 1990 June 29. This striking image showed 
an X-ray-bright sunlit half-circle on one side, and a much 
dimmer but not completely dark side outlined by a brighter 
surrounding diffuse X-ray background. Several origins for the 
dark-side emission were considered, but the authors' favored 
explanation was continuum emission arising from solar wind 
electrons sweeping around to the unlit side and impacting 
on the lunar surface, producing thick-target bremsstrahlung. 
Given the very limited energy resolution of the PSPC, how- 
ever, emission from multiple lines could not be ruled out. 

A significant problem with the bremsstrahlung model was 
explaining how electrons from the general direction of the Sun 
could produce events on the opposite side of the Moon, with a 
spatial distribution which was "consistent with the telescope- 
vignetted signal of a constant extended source." An elegant al- 
ternative explanation would be a source of X-ray emission be- 
tween the Earth and the Moon, but at the time, no such source 
could be envisioned. If this source were also time-variable, it 
would account for the Long Term Enhancements (LTEs) seen 
by ROSAT . These occasional increases in the counting rate of 
the PSPC are vignetted in the same way as sky-b ackground 
X-rays, indicating an external origin ( Snowden et al. 1995). 
LTEs are distinct from the particle-induced background, and 
are uncorrelated with the spacecraft's orientati on o r posi- 
tion (geomagnetic latitude, etc.), although Frevberg (1994) 
noted that LTEs appeared to be related, by a then unknown 
mechanism, to geomagnetic storms and solar wind variations. 
The fi nal ROSAT All-Sky Survey ( RASS) diffuse background 
maps (Snowden et al. 1995, 1997) removed the LTEs, so far 
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as possible, by comparing multiple observations of the same 
part of the sky, but any constant or slowly varying (r > 1 
week) emission arising from whatever was causing the LTEs 
would remain. 

A conceptual breakthroug h came with the R OSAT obser- 
vation of c omet Hyakutake tLisse et al.lll996l) and the sug- 
gestion by ICravensI i 19971) that charge transfer (CT) be- 
tween the solar wind and neutral gas from the comet gave 
rise to the observed X-ray emission. In solar-wind charge 
transfer, a highly charged ion in the wind (usually oxygen 
or carbon) collides with neutral gas (mostly water vapor 
in the case of comets) and an electron is transferred from 
the neutral species into an excited energy level of the wind 
ion, which then decays and emits an X-ray. This hypoth- 
esis has been proven by subsequent obser vations of comets 
such as C/LINEAR 1999 S 4 by Chandra dLisse et alJ pOOl fc 
and Hyakutake by EUVE ( Krasnopolskv & Mumma 2001) 
(see also the review by ICravensI J2002I) ), and is supp orted 
by in c reasingly detailed spectral mod els dKharchenko et alJ 
2003; Kharchenko & Dalgarno 2000). A more extensive 
history of the evolution of the solar- wind CT conce pt can 
be found in ICra vens. R obertson. & Snowdenl ( 12001b and 
iRobertson & Cravens! (I2003all. 

Citing the cometary emission model, ICoxl JT998h pointed 
out that CT must occur throughout the heliosphere as the 
solar wind interact s with atomic H and He within the solar 
system. Frevberg ( 1998) likewise presented ROSAT High- 
Resolution Imager data that provided some evidence for a cor- 
relation between increases in the apparent intensity of comet 
Hyakutake and in the detector background; he further sug- 
gested that this could be caused by charge transfer of the so- 
lar wind with the Earth's atmosphere. A rough broad-band 
quantitative analysis by ICravensI (120001) predicted that he- 
liospheric emission, along with CT between the solar wind 
and neutral H in the Earth's tenuous outer atmosphere (geo- 
corona), accounts for up to half of the observed soft X-ray 
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TABLE 1. 
Observation Information 









Exposure 


Chandra 


ObsID 


Date 


CCDs 


(s) 


Time 


2469 


2U01 Jul 26 


123, S23 


2930 


1 125000 /0-1 12503000 


2487 


2001 Jul 26 


123, S23 


2982 


112503320-112506302 


2488 


2001 Jul 26 


123, S23 


2747 


112507858-112510605 


2489 


2001 Jul 26 


123, S23 


2998 


112510900-112513898 


2490 


2001 Jul 26 


123, S23 


2830 


112515450-112518280 


2493 


2001 Jul 26 


123, S23 


2993 


112518500-112521493 


2468 


2001 Sep 22 


123, S123 


3157 


117529483-117532640 


3368 


2001 Sep 22 


123, S123 


2223 


117532850-117535073 


3370 


2001 Sep 22 


123, S123 


4772 


117536678-117541450 


3371 


2001 Sep 22 


123, S123 


3998 


117541880-117545878 


Note. 


— Chandra time corresponds to the beginning of 1998. Obser- 



vations run from July 26 02:01:10-07:58:13 UT, and September 22 07:04:43- 
11:37:58 UT. 



background (SXRB). Intriguingly, results from the Wisconsin 
Soft X -Ray Background sky survey (McCammon & Sanders 
1990) and RASS observations l lSnowden et all 1995ft indicate 
that roughly half of the 1/4-keV background comes from a 
"local hot plasma." Cravens. Rob ertson. & Snowdenl (120011) 
also modeled how variations in solar-wind density and speed 
should affect heliospheric and geocoronal CT emission ob- 
served at Earth, and found strong correlations between the 
measured solar-wind proton flux and temporal variations in 
the ROSAT counting rate. 

In this paper we present definitive spectral evidence for 
geocoronal CT X-ray emission, obtained in Chandra obser- 
vations of the Moon. Data analysis is discussed in |2] and 
results are presented in S|3] As we show in ^ model pre- 
dictions of geocoronal CT agree very well with the observed 
Chandra spectra. In *J5]we estimate the level of heliospheric 
CT emission, discuss the overall contribution of CT emission 
to the SXRB, and assess the observational prospects for im- 
proving our understanding of this subtle but ubiquitous souce 
of X-rays. 

2. THE DATA 

The Moon was observed with the Chandra Advanced CCD 
for Imaging Spectroscopy (ACIS) in two series of calibration 
observations on 2001 July 26 and September 22 totaling 17.5 
and 14 ksec, respectively (see Table 0. The intention was 
to determine the intrinsic ACIS detector background by us- 
ing the Moon to block all cosmic X-ray emission. Four of 
the ACIS CCDs were used in July (12 and 13 from the ACIS- 
Imaging array, and S2 and S3 from the ACIS-Spectroscopy 
array), and the S 1 chip was added in September. Two of the 
chips, SI and S3, are back-illuminated (BI) and have better 
quantum efficiency at low energies than the front-illuminated 
(FI) chips, 12, 13, and S2. As can be seen in Fig. ^ how- 
ever, the BI chips have higher intrinsic background than the 
FI chips, and also poorer energy resolution. Telemetry lim- 
its prevented the operation of more CCDs when using ACIS 
Very Faint mode, which was desired because of its particle- 
background rejection utility ( Vikhlinin 2001). 

The ACIS detector background was also calibrated in 
an alternative manner usin g Eve nt Histogram Mode (EHM; 
iBiller. Plucinskv. & Edgar] ( 120021) '). The July Moon and EHM 



spectra from the S3 ship were compared by Mark evitch et alJ 
(2002) and showed good agreement, although there was a no- 
ticeable but statistically marginal excess near 600 eV in the 
Moon data. 

The dark-Moon vs EHM comparison strongly supports the 
assumption that the high-energy particle background inside 
the detector housing where EHM data are collected is the 
same as in the focal position. With that in mind, new cali- 
bration measurements were made on 2002 September 3 with 
ACIS operating with its standard imaging setup in a "stowed" 
position, where it was both shielded from the sky and removed 
from the radioactive calibration source in its normal off-duty 
position. dMarkevitohll2002l) . These data (ObsID 62850, 53 
ks) provide the best available calibration of the intrinsic de- 
tector background and are used in the analysis that follows. 

2.1. Data Preparation 

All data were processed to level 1 using Chandra Interac- 
tive Analysis of Observations (CIAO) software, Pipeline re- 
lease 6.3.1, with bad-pixel filtering. Start and stop times for 
each observation were chosen to exclude spacecraft maneu- 
vers. Apart from the inclusion of CTI corrections (see below) 
and a more aggressive exclusion of any possibly questionable 
data, our data processing is essent ially the same as that de- 
scribed by Mar kevitch et alJ (|2002), who limited their analy- 
sis to the July S3 data. Here we use data from all chips during 
both the July and September observations, and include data 
from periods of pa rtial dark-Moon coverage by using spatial 
filtering (see 32.2I >. 

Although ACIS has thinly aluminized filters to limit optical 
contamination, the sunlit side of the Moon was so bright that 
an excess bias signal was sometimes produced in the CCDs, 
particularly in the 12 and 13 chips that imaged that region dur- 
ing July. As described in Markevit ch et alJ ( 120021) . a bias cor- 
rection to each event's pulse height amplitude was calculated 
by averaging the 16 lowest-signal pixels of the 5x5-pixel 
Very Faint (VF) mode event island. ObsID 2469 suffered by 
far the most optical contamination, so that all data from the 12 
chip during that observation had to be discarded. The 13 chip 
also had significant contamination, but it was s mall enough 
to be largely corrected. As explained in 32.31 however, 13 
data from that observation were also excluded as a precau- 
tion. Two other July observations required exclusion of some 
12 data because of optical leaks (930 s in ObsID 2490 and 
1 140 s in Obsid 2493), but in all the remaining data the typ- 
ical energy correction was no more than a few eV, which is 
insignificant for our purposes. 

To improve the effective energy resolution, we applied stan- 
dard charge transfer inefficiency (CTI) corrections, as imple- 
mented in the CIAO tool acis_process_event s, to data 
from the FI chips. CTI is much less of a problem in the 
BI chips, SI and S3, and no corrections were made to those 
chips' data . Finally, V F-mode filtering was applied to all the 
data ( I Vikhlininl 1200 lh to reduce the particle-induced detec- 
tor background. The "ACIS stowed" background data were 
treated in the same way, except that no optical-contamination 
corrections were required. 

2.2. Spatial Extractions 

As seen in Fig. ^ Chandra pointed at a fixed location on 
the sky during each observation as the Moon drifted across 
the field of view. Using ephemerides for the Moon and Chan- 
dra's orbit, we calculated the apparent position and size of 
the Moon in 1 -minute intervals and extracted data from its 
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dark side, as well as from the bright side and from the un- 
obscured cosmic X-ray background (CXB) within the field of 
view for comparison. Because the Moon moved by up to 16" 
per minute, and to avoid X-ray "contamination" of data within 
each extraction region (particularly spillover of bright-side or 
CXB photons into the dark side) we used generous buffers of 
90" from the terminator and limb for the dark-side extraction, 
30" from the terminator and 60" from the limb for bright-side 
data, and 60" from the limb for the CXB data. Chandra has 



July 2001 



Obsid 2469 


Obsid 2487 




(^^^/ 


Obsid 2488 


Obsid 2489 


'' 'dlk^A 




Obsid 2490 


Obsid 2493 






September 2001 


Obsid 2468 


Obsid 3368 






Obsid 3370 


Obsid 3371 



FIG. 1. — Moon motion across the ACIS chips during July and September 
observations. Green denotes Moon position at the start of the observation 
and red at the end. The illuminated portion of the Moon is the crescent on 
the right. Dark-side gibbous extraction regions provide a 90" buffer along 
the Moon limb and terminator. SI and S3 chips are darker because of their 
higher background. 



a very tight point spread function, with an on-axis encircled 
energy fraction of nearly 99% at 500 eV within 10"; although 
off-axis imaging is involved here, estimated X-ray contami- 
nation is less than ~ 2% within our chosen extraction regions. 

Data from all ObsIDs were analyzed in several energy 
bands to look for discrete sources but none were found, and 
lightcurves for each observation behave as would be expected 
for uniform emission within each extraction region. Effective 
exposure times (as if one CCD were fully exposed) were com- 
puted for each ObsID/chip/extraction combination by com- 
puting extraction areas for each 1 -minute interval (accounting 
for spacecraft dither, which affects area calculations near the 
chip edges) and summing the area x time products. Results are 
listed in Table|2] 

Because the detector background is not perfectly uniform 
across each chip, background data were projected onto the 
sky and extracted using the same regions as for the dark- 
Moon data. Exposure-weighted and epoch-appropriate de- 
tector response functions (RMFs and ARFs) were then cre- 
ated using standard CIAO threads, including the corrarf 
routine, which applies the ACISABS model to account for 
contaminant build-up on ACIS. The detector background rate 
varies slightly on timescales of months, so we renormalized 
the background data to match the corresponding observational 
data in the energy range 9.2-12.2 keV, where the detected sig- 
nal is entirely from intrinsic background. The required adjust- 
ments were only a few percent. 

2.3. Spectra 

As described by M arkevitch et alJ J2002I) . the BI chips, and 
very rarely the FI chips, often experience "soft" background 
flares because of their higher sensitivity to low-energy parti- 
cles. A relatively bright flare was found in ObsID 3370, and 
400 seconds of data were removed. Weaker flares are more 
common and we judged it better to model and subtract their 
small effects rather than exclude large intervals of data. Soft 
flares have a consistent spectral shape (a powerlaw with high- 
energy cutoff) and their intensity at all energies can be deter- 
mined by integrating the excess signal (above the "stowed" 
background) in the energy range 2.5-7 keV where the relative 
excess is most significant. We find that spectra from ObsIDs 
2468 and 3368 have minor soft-flare components, and have 
accounted for them in the results presented later. An essen- 
tially negligible soft-flare excess is also seen and accounted 
for when all the July S3 data are combined. 

One last complication is the effect of optical contamination 
on the energy calibration. Although raw energy offsets were 
removed from the data during the bias-correction process, 
more subtle effects remained, mostly related to charge trans- 
fer inefficiency. Optical-leak events partially fill the charge 
traps in the CCDs, thus reducing the net CTI. When standard 
CTI corrections were applied to the data to improve the en- 
ergy resolution of the FI chips, this overcorrected and pushed 
energies too high for data affected by the optical leak. The 
effect varies between and within chips based on optical expo- 
sure, but we can place an upper limit on it by examining the 
bright-Moon data, which are most affected. Fig.|2]shows the 
spectrum of the combined July bright-Moon data, which come 
from the 12 and 13 chips. The /T-shell fluorescence lines of O, 
Mg, Al, and Si are easily identified, and we find a positive 
offset of 50 ± 5 eV from their true values of 525, 1254, 1487, 
and 1740 eV, respectively. Bright-side data from ObsID 2469 
13 were excluded because they showed an offset of roughly 80 
eV, with a distorted shape for the O peak; to be conservative, 
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TABLE 2. 

Effective Exposure Times Per Chip" 







Dark-Side Region 






Bright-Side Region 




WDSl 




13 


SI 


S2 




T9 
L— 


13 


SI 


S2 




2469 


o b 


C 




2930 


2924 


o b 


C 










2487 


1702 


1759 




2812 


2982 


722 


704 










2488 


1204 


1233 




2714 


2747 


979 


965 




o 





2489 


1491 


1375 




2978 


2977 


984 


1062 










2490 


1071 d 


1306 




2811 


2714 


984 


1062 










2493 


1456 e 


1637 




2862 


2572 


238 


847 










Total 


6924 


7310 




17107 


16916 


3503 


4560 










2468 


3157 


3150 


1026 


2534 


2958 








1462 


177 


6 


3368 


2223 


2223 


162 


1364 


2031 








1680 


275 





3370 


4772 


4770 


1281 f 


3479 


3611 f 








2157 f 


439 


20' 


3371 


3995 


3995 


265 


1631 


1803 








3024 


1114 


950 


Total 


14147 


14138 


2734 


9008 


10403 








8323 


2005 


976 



" In units of s. * Severe optical leak. All data excluded. ' Unreliable energies and possible event 
loss from optical leak. All data excluded. d Optical leak. Time range 1 12516560-1 12517490 (930 
s) excluded. * Optical leak. Time range 1 12520070-1 12521210 (1 140 s) excluded. f Background 
flare in BI chips. Time range 1 17538500-1 17538900 (400 s) excluded. 



and the agreement of astrophysical line positions in the F I and 
S3 spectra with each other and theoretical models (see $4.11 . 
the dark-side energy errors indeed appear to be negligible. 

3. RESULTS 
3.1. Dark-Side Spectra 

X-ray spectra were created from the event files using CIAO 
dmextract. Data from the three FI chips (12, 13, and S2), 
which have lower QE than the BI chips below 1 keV, were 
always combined in order to improve statistics. Scaled back- 
ground spectra, with soft-flare corrections as needed, were 
created for each observed dark-side spectrum and subtracted 
to reveal any excess X-ray emission. 

Summing all the July data apart from ObsID 2469 reveals 
no significant excesses in either the S3 or combined FI spectra 
(see top half of Fig.|5Jl. The S3 spectrum from ObsID 2469, 
however, has a noticeable emission feature (more than 3a) at 
~ 600 eV. The corresponding FI spectrum for ObsID 2469 
has too few counts to be used for corroboration as it includes 
only data from the S2 chip (because of the optical leaks in 12 
and 13). 

Much stronger evidence for excess emission near 600 eV 
appears in the September dark-Moon spectra, in all chips (bot- 
tom half of Fig. [3}- It is immediately obvious that this can not 
be particle- or photon-induced O fluorescence, which would 
occur in a single line at 525 eV, nor is it electron-impact con- 
tinuum as posited by Schmitt et alj jl991l> . 

To assess the significance of any excesses, we selected three 
energy ranges for statistical study with the a priori assumption 
that solar-wind charge transfer is the source of the emission 
($0. Each range (311-511 eV, 51 1-71 1 eV, and 716-886 eV) 
was chosen to extend ~ 50 eV below a nd a bove the strongest 
CT lines expected within each range ( 34. 1> . For comparison, 
the 886-986 eV band was also studied (in which we might 
hope to see Ne IX Ka at 905-922 eV), along with four 200- 



we excluded the corresponding dark-side data from further 
consideration as well. As noted before, 12 data from ObsID 
2469 were already excluded because of their much larger op- 
tical contamination. 

Energy offsets in the dark-side and CXB data should be 
much smaller, particularly for the S chips, which did not im- 
age the bright side of the Moon. Judging from the positions 
of weak fluorescence lines present in the detector background, 
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FIG. 2. — Spectrum of the bright side of the Moon, combining 12 and 13 
data from all July observations except ObsID 2469, binned by 2 PI channels 
(29.2 eV). Dotted curve is detector background. Fluorescence K-shell lines 
from O, Mg, Al, and Si are shifted up by 50 eV from their true values because 
of residual errors when correcting for detector sensitivity to optical photons 
(see text). Optical contamination effects likewise cause slight mismatches in 
energies of intrinsic detector features such as the Au-M complex (2.2 keV), 
Ni-K (7.5 keV), and Au-La (9.7 keV). There are ~1300 counts in the O-K 
line. 
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Sept 2001, FI chips 
- 1 — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — 

Total observed 

Background 




2500 




1500 



2000 



1000 1500 
Energy (eV) 



2000 



2500 




2500 



FIG. 3. — Observed and background-subtracted spectra from the dark side of the Moon, with 2-channel (29.2 eV) binning. In the July S3 data (top four 
panels), no excess emission is obvious except in ObsID 2469. Features around 1750 eV (Si-K) and 2150 eV (Au-M) are from particle-induced fluorescence of 
the detector assembly. Bottom four panels show S3 and combined FI spectra from the three September ObsIDs with the strongest emission excesses. Oxygen 
emission from charge transfer is clearly seen in both spectra, and energy resolution in the FI chips is sufficient that O VIII Lya is largely resolved from O VII 
Ka. High-n O VIII Lyman lines are also apparent in the FI spectrum, along with what is likely Mg XI Ka around 1340 eV. 



eV-wide bands from 1000 to 1800 eV. The most important 
results are shown in Table[3] with excesses of more than 2.5a 
shown in bold. 

Emission in the 511-711 eV range has an excess of more 
than 4<t in three of the four S3 spectra from September, and 
2.3cr in the other (ObsID 3370). The same pattern (most sig- 
nificant excess in ObsID 2468, least in ObsID 3370) holds for 
the combined FI spectra. When spectra from the three ObsIDs 
showing the largest excesses (2468, 3368, and 3371, hence- 
forth referred to as the "bright-September" ObsIDs) were 
summed, the feature significance was more than 6a in both 
the S3 and FI spectra (see Table The ratio of net count- 
ing rates for S3 and the FI chips (4.3 ±0.9) also matches well 
with the ratio of those chips' effective areas in that energy 



range, consistent with this being an X-ray signal from the sky. 
Results from the SI chip were consistent with those from S3, 
but with lower significance because of the Si's much shorter 
dark-Moon exposure times and somewhat higher background; 
we do not discuss S 1 dark-side results further. 

The 716-886 eV range, which we expect to contain O VIII 
Ly/3 and O VIII LyS emission, also showed significant ex- 
cesses in the S3 and combined FI spectra for the bright- 
September ObsIDs (~ 2.5cr and ~ 4a, respectively) with ex- 
cesses in individual ObsIDs roughly following the time pat- 
tern seen in the 51 1-71 1-eV band. The S3 spectrum for Ob- 
sID 3371 stands out with a 4a excess. The same ObsID also 
has 2a excesses in the 886-986 eV range in both the S3 and FI 
spectra. This energy range, like the 716-886 eV range, con- 
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TABLE 3. 

Net Emission Within Selected Energy Bands 



311-511 eV 511-716 eV 716-886 eV 886-986 eV 



Counting Signif. Counting Signif. Counting Signif. Counting Signif. 



Chips 


ObsIDs 


Rate* 




Rate a 


(v) 


Rate" 


(<r) 


Rate" 


W) 


FIs" 


2469 


2± 19 


0.1 


25 ±20 


1.2 


6± 16 


0.4 


25 ±16 


1.5 


FIs 


all July c 


-1±6 


-0.2 


3±6 


0.5 


7±5 


1.5 


3±4 


0.7 


FIs 


2468 


2±12 


0.1 


65 ±14 


4.5 


32 ±12 


2.8 


-1±7 


-0.1 


FIs 


3368 


-11 ± 13 


-0.9 


76 ±18 


4.2 


22 ± 13 


1.7 


18±11 


1.7 


FIs 


3370 


-5 ±9 


-0.5 


12±9 


1.3 


18±9 


2.1 


2±6 


0.3 


FIs 


3371 


1 ± 11 


0.1 


26 ± 12 


2.2 


22 ± 10 


2.1 


15 ± 8 


1.9 


FIs 


3 bright d 


-2 ±7 


-0.3 


52±8 


6.3 


26±7 


3.9 


10±5 


2.1 


S3 


2469 


-36 ±34 


-1.1 


126 ±39 


3.2 


29 ±30 


1.0 


-19 ± 17 


-1.1 


S3 


quiet July c 


30±18 


1.7 


22 ± 15 


1.4 


13 ± 13 


1.0 


6±10 


0.7 


S3 


all July c 


22 ±16 


1.4 


43 ±14 


3.0 


18 ± 12 


1.5 


4±8 


0.4 


S3 


2468 


51 ±31 


1.7 


218 ±33 


6.6 


27 ±23 


1.2 


9± 16 


0.5 


S3 


3368 


54 ±54 


1.0 


256 ± 61 


4.2 


18±40 


0.4 


11 ±30 


0.4 


S3 


3370 


58 ±37 


1.5 


77 ±34 


2.3 


-21 ±24 


-0.9 


-15 ±17 


-0.9 


S3 


3371 


26 ±49 


0.5 


247 ± 60 


4.1 


214 ±54 


4.0 


72 ±35 


2.0 


S3 


3 bright d 


63 ±28 


2.2 


228 ±31 


7.3 


58 ±23 


2.6 


18±16 


1.2 



NOTE. — Rate excesses of more than 2.5cr are shown in bold. "FIs" means the 12, 13, and S2 chips in combination. 
" Units of 10~ 6 cts s~' arcmin~ 2 . b Chip S2 only; data from 12 and 13 are unusable. c ObsIDs 2469, 2487, 2488, 2489, 
2490, and 2493. d The "bright-September" ObsIDs: 2468, 3368, and 3371. e All July ObsIDs except 2469. 



tains CT emission lines (Ne IX Ka at 905-922 eV) from an 
ion which is only abundant when the s olar wind is especially 
highly ionized. As we will discuss in 34.31 there is evidence 
for such a situation during ObsID 3371. 

Above 1000 eV, no significant excesses were seen for any 
of the ObsID combinations listed in Table|3]with the possible 
exception of 1200-1400 eV, which had a 2.9er excess in the 
bright-September FI spectrum (and a l.Ocr excess in the S3 
spectrum). Again, 3371 was the individual ObsID recording 
the largest excesses in the FI (1.8<r) and S3 (1.3a) spectra. 
Although the evidence is not compelling, we believe that the 
observed excess probably represents a detection of He-like 
MgXlKa(- 1340 eV). 

In the 3 1 1-5 1 1 eV range where C VI Lyman emission might 
be detectable, only a 2.2cr excess appears in the S3 data. Over 
the full range of O emission (51 1-886 eV), which is relevant 
to the discussion in 50] the bright-September S3 and FI spec- 
tra both have an excess of 7.4<j, with net rates of 287 ± 39 and 
78 ±11 x 10~ 6 cts s -1 arcmin -2 , respectively. The summed 
July S3 data have a rate of 62 ± 19 x 10~ 6 cts s" 1 arcmin -2 
between 511 and 886 eV. If ObsID 2469 is excluded because 
of its obviously stronger O emission, the rest of the July S3 
data have a statistically insignificant 1 .7 a excess with a rate 
of 34 ±20 x 10" 6 cts s -1 arcmin -2 . 

(As an aside, we note that the S3 Moon spectrum from 
July, combining data from all six ObsID s, was used as a mea - 
sure of the detector background by Markevi tch et alJ |2002). 
Even with the ObsID-2469 excess, the oxygen emission in 
that background is much less than the sky-background emis- 
sion discussed in that paper, and so the authors' results are not 
significantly affected.) 



3.2. Comparison with CXB 

If the dark-side emission arises between the Earth and 
Moon, then it must also be present at the same level on the 
bright side and in the cosmic X-ray background beyond the 
Moon's limb. Unfortunately, the bright side was observed us- 
ing the 12 and 13 chips and only in July, when the dark-side 
emission was barely detectable even with the more sensitive 
S3 chip. Such a weak signal would in any case be swamped 
by the bright-side fluorescence X-rays. 

Spectra of the CXB were obtained only in September, pri- 
marily by the S 1 chip (see Table[2} which is back-illuminated 
like S3 and has similar quantum efficiency. In Fig. [4] which 
shows CXB S 1 spectra from all four September observations, 
it is apparent that the CXB is much brighter than the dark- 
Moon signal. In fact, this is one of the brightest regions of 
the sky, with a complex spatial structure; see Table 0] which 
lists the centroid of each CXB extraction region and the corre- 
sponding R45 (3/4-keV-band) RASS rate. 3 The typical SXRB 
recorded in the RASS is roughly one-quarter as bright, com- 
parable to the bright-September dark-Moon brightness. 

One can also see that emission around 600 eV is strongest 
for ObsID 2468, which is also the ObsID that shows the most 
significant dark-side X-ray emission. While this correlation 

3 Background maps are available online at 

http://www.xray.mpe.mpg.de/rosat/survey/sxrb/12/ass.html 
An X-ray Background Tool is also available at 
http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pl Given the 
limited resolution of the PSPC, data are usually divided into three energy 
bands: R12 (a.k.a. the 1/4-keV band, effectively defined on the high-energy 
end by the C absorption edge at 0.284 keV), R45 (3/4-keV band; roughly 
0.4-1.0 keV), and R67 (1.5-keV band; roughly 1.0-2.0 keV). 
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TABLE 4. 

ROSAT All-Sky Survey Background Rates for CXB 
Regions 









R45 Rate 


ObsID 


Centroid RA 


Centroid Dec 


(10 -6 cts s -1 arcmin~ 2 ) 


2468 


15*41 m 44 I 


-37°00'30" 


543.7 ± 103.2 


3368 


15*41'"29 ! 


-37° 15' 10" 


483.3 ±100.3 


3370 


15''41'"18 s 


-37°29'10" 


549.0 ± 107.1 


3371 


15 A 41"'30 S 


-37°41'20" 


468.7 ±100.5 


Note. 


— R45 band 


is defined as 


PI channels 52-90, cor- 



responding to approximately 0.4-1.0 keV. ROSAT rates were 
found using the HEASARC X-Ray Background Tool version 2.1 
(http://heasarc.gsfc.nasa.gov/Toolsl with a cone radius of 0.1 de- 
grees. 



is suggestive, direct comparisons among the CXB spectra are 
not possible because they were taken from slightly different 
regions of the sky (because of the Moon's motion), nor can 
the CXB data be adequately normalized using ROSAT All- 
Sky Survey background rates because statistical uncertainties 
are too large. 

4. INTERPRETATION: GEOCORONAL CHARGE TRANSFER 

The primary result of our analysis is that highly-significant 
time-variable emission is seen looking toward the dark side of 
the Moon at energies between 500 and 900 eV. As we discuss 
in this section, the observed spectrum, intensity, and temporal 
behavior can all be explained by charge transfer between the 
solar wind and the Earth's outer atmosphere. 

Charge transfer is the radiationless collisional transfer of 
one or more electrons from a neutral atom or molecule to an 
ion. When the recipient ion is highly charged, it is left in a 
high-n excited state which then decays via single or sequential 
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FIG. 4. — Background-subtracted CXB spectra from the SI chip for all 
four September observations, with 3-channel binning (43.8 eV) for clarity. 
Spectra are not directly comparable because they were collected from slightly 
different regions of the sky (see Table Bl. all of which are quite bright. The 
bright-September dark-side S3 spectrum is shown for comparison. 



TABLE 5. 

Slow Solar Wind Ion Abundances and Cross 
Sections 





Abund. 


acr with H 


aci with He 


Ion 


rel. to O 


(10~ 15 cm 2 ) 


(10~ 15 cm 2 ) 


C VI 


0.21 


1.03 


1.3 


C VII 


0.32 


4.16 


1.5 


NVII 


0.06 


3.71 


1.7 


N VIII 


0.006 


5.67 


2.0 


OVII 


0.73 


3.67 


1.1 


O VIII 


0.20 


3.4 


1.8 


OIX 


0.07 


5.65 


2.8 


NelX 


0.084 


3.7 


1.5 


NeX 


0.004 


5.2 


2.4 



NOTE. — Relative element abundances: C/O = 0.67, N/O 
= 0.0785, Ne/O = 0.088. 



radiative transitions. (The ion may also autoionize if multiple 
electrons are transferred from the neutral.) In geocoronal CT, 
the neutral gas is atomic H in the Earth's outer atmosphere 
extending tens of thousands of km into space, the X-ray emit- 
ting ions are heavy elements such as C, O, and Ne in the solar 
wind, and the collisions take place outside the magnetosphere, 
into which the wind particles can not penetrate. At X-ray en- 
ergies, most of the emission comes from hydrogenic and He- 
like C and O ions because of their relatively high abundance. 

4. 1 . Model Spectra 

The equation for the CT emissivity of line / from ion i can 
be written as 

ea = v c n n niy a <Ji phot s _I cm" 3 , (1) 

where v c is the collision velocity (effectively the solar wind 
velocity), n„ is the neutral species density, «,■ is the relevant 
ion density, yu is the net line emission yield per CT-excited 
ion, and er, is the total CT cross se ction for ion 

As described in the review by ISmith et alJ J2003), solar 
wind velocity and ionization level are closely correlated, and 
can be used to divide the solar wind into two main types: a 
"slow" highly ionized wind (v f ~ 400 km s" 1 ), and a "fast" 
less ionized wind (v c ~ 700 km s" 1 ). During solar minimum, 
the slow wind is found near the equatorial plane (solar lati- 
tudes between roughly -20° and +20°), while the fast wind 
dominates at higher latitudes. During solar maximum, the 
slow wind extends to higher latitudes, but there is significant 
mixing of the two components. Avera ge slow-wind ion abun- 
dances relative to oxygen, taken from Schw adron & Cravens! 
(2000), are listed in Table[5] along with cross sections for CT 
with neutral H and He. 

A great deal of physics is contained within yu<Ji, as it in- 
cludes the initial quantum-sublevel population distribution of 
the ion immediately following electron capture, and then the 
branching ratios from all those levels during the subsequent 
radiative cascades. Cross sections for CT of highly charged 
C, N, O, and Ne ions with ato mic H and associated r adiative 
branching ratios are taken from Harel, Jouin, & Pons (1998), 
iGrjsgnwood^ydJ l|2001l). |Rigazio. Kharchenko. & Dalgarno 
(2002), Johns on et alJ d2000t). and related references in 
Kharchenko & Dalgarno (2000, 2001). Total CT cross sec- 
tions for all ions are a few xl0~ 15 cm" 2 , with uncertainties 
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TABLE 6. 
Major Solar CT Emission Lines 



TABLE 7. 

Fit Results versus Model Line Brightness 



Line 


Energy 
(eV) 


Line Yield 


C V Ka a 


299, 304, 308 


0.899 


CVILya 


368 


0.650 


N VI Ka a 


420,426, 431 


0.872 


C VI Ly/3 


436 


0.108 


C VI Ly7 


459 


0.165 


VII Ka a 


561, 569, 574 


0.865 


VIII Lya 


654 


0.707 


VII K/3 


666 


0.121 


VIII Ly/3 


775 


0.091 


VIII Ly7 


817 


0.033 


O VIII Ly5 


836 


0.103 


O VIII Lye 


847 


0.030 


Ne IX K« a 


905, 916, 922 


0.887 



"The He-like Ka complex consists of three lines: forbidden 



° 2 'So), intercombination (\s2p i P\,i 



Is 1 'So), 



and resonance (\s2p 'P\ — > Is? 'So). For charge transfer, the 
forbidden line (lowest energy) dominates. 



of typically 30%, and are fairly constant as a function of col- 
lision velocity near 400 km s" 1 . Cross sections for electron 
capture into individual sublevels have larger errors, which are 
the major contributors to uncertainties in the line yields, par- 
ticularly for H-like ions. As can be seen from the line yields 
listed in Table|6j emission from He-like ions is predominantly 
(~90%) in the form of Ka (n = 2 — » 1) photons, while H-like 
emission is split more evenly between Lya and the higher-n 
transitions (e.g., Ly7 and Ly£). The unusual strength of the 
high-n Lyman lines is a unique signature of CT which can not 
be reproduced by thermal plasmas. 

For comparison with wo rk by 

Cravens, Robertson, & Snowden and 
Robertson & Cravens (2003a), we calculate the value of 
a = J2u( n i/ n p)yn (J iEii, where n p is the solar-wind proton 
density, En is the line energy, and the sum is over all CT 
lines from 9 5 eV (the lower limit o f ROSAT 's energy range) 
to 1000 eV. iRobertson & Cravens! J2003al) estimate that a 
equals 6 x 10" 16 eV cm 2 for CT with H and 3 x 10" 16 eV 
cm 2 for He. We derive 9.82 x 10" 16 and 4.56 x 10" 16 eV 
cm 2 , respectively, or 8.13 x 10" 16 and 3.92 x 10" 16 eV cm 2 
for energies above 180 eV, where the ROSAT PSPC effective 
area becomes appreciable (see Fig. Whatever its value, 
a global parameter such as a is insufficient for the work 
described here; when analyzing data with at least moderate 
energy resolution it is necessary to create a line emission 
model. 

Our resulting model spectrum for geocoronal CT is shown 
in Fig. [5] That model was then used to simultaneously fit 
the bright-September S3 and FI data between 250 and 5000 
eV using their associated ACIS response functions (see 32. 2> . 
Given the limited statistics (~ 700 counts attributable to X- 
ray emission in the S3 data, and fewer in the FI spectra), 
we grouped the model emission into four "lines" at 440, 570, 
660, and 810 eV; the 440-eV and 810-eV lines are modeled 
as finite-width Gaussians since they represent several lines 
spread over a 100-eV and 60-eV range, respectively. 



Energy 
(eV) 


Range a 
(eV) 


Major 
Lines 


B/it 


R c 

D mode! 


B fj, 

"model 


440 


367^170 


C Lya-5 


32 ( 0-76) 


5.58 


6 


570 


561-574 


OKq 


85 (55-115) 


2.38 


36 


660 


654-666 


O Lya, K/3 


38 (23-52) 


0.96 


40 


810 


775-847 


O Ly/3-e 


15 (9-22) 


0.23 


66 


1340 




Mg Kq 


4(1-7) 







a Range of energies for grouped CT-model lines. 

b Fit results for the 3 bright-September ObsIDs, in units of 10~ 



phot 



c Model line brightnesses are calculated using the average slow- 
solar-wind parameters listed in the text, in units of 10~ 8 phot s~ cnf 2 
arcmin~ 2 . Detailed model predictions for the Mg XI Ka line at 1340 
eV are not available, but its normal brightness is negligibly small. 



Because of the poor statistics, the only free parameter is 
the normalization. As can be seen in Fig. |6j overall agree- 
ment between the shape of the model spectrum and observa- 
tion is good, particularly for the O lines which are the most 
prominent features. The fit shows that carbon emission is 
overpredicted, however, which may be because of atypical ion 
abundances, errors in CT cross sections and line yields, or un- 
certainties in the QE of the ACIS detectors near the extreme 
limits of their effective energy ranges. Better agreement is ob- 
tained when relative ion abundances are adjusted; in the ad- 
justed fit, we reduce the abundances (relative to He-like O VII) 
of H-like C VI and O VIII by factors of 6 and 2, respectively. 
When a line is added to fit the putative Mg XI Ka feature at 
1340 eV, the F-test significance is 99.6%. 

In fitting the FI and S3 spectra simultaneously, we are im- 
plicitly assuming that they recorded the same emission. This 
is not strictly true because of differing exposures for each 
chip/ObsID combination, and in fact, if the S3 and FI spec- 
tra are fit separately, the S3 brightness comes out somewhat 
higher than for the FI data. Although within statistical un- 



500 




300 



400 600 
Energy (eV) 

FIG. 5. — Model geocoronal X-ray spectrum and effective areas for Chan- 
dra and ROSAT . Brightness is for an observation through the magnetosphere 
flank assuming average slow-solar-wind parameters, plotted with FWHM en- 
ergy resolution of 6 eV. The cluster of lines below 150 eV comprises mostly 
n = 4, 5 — > 2 transitions in O VI and O VII. 
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FIG. 6. — Background-subtracted September spectra (3 brightest ObsIDs) and charge transfer model fits using four composite lines. The same model, including 
normalization, is used to fit the S3 and FI spectra. The nominal fit uses the average (baseline) solar-wind parameters discussed in the text. In the adjusted fit, the 
relative C VI vs O VII Ka emission is 1/6 of the nominal case, and O VIII emission is reduced by half. Results for the adjusted fit are listed in Tablel71 



certainties, this difference is consistent with the sense of sus- 
pected errors in the quantum efficiency of the ACIS chips. 

Line brightnesses from the adjusted-abundances fit are 
listed in Table0 note that O VII Ka and O VIII Lya are essen- 
tially unresolved in the S3 data so their combined brightness 
is more reliable than the individual values. We also list the 
corresponding model pred ictio ns for an average solar wind 
(which will be derived in jj4.2t . which are more than an or- 
der ofmagnitude smaller than observed. As we will discuss 
in 34.31 based on available solar-wind data, there is good rea- 
son to expect a much higher-than-average geocoronal emis- 
sion level during the September observations, as well as rela- 
tively weak C VI emission. 

We lastly note that the bright-September S3 spectrum is 
remarkably similar in shape (and normalization) to the dif- 
fuse X -ray background spectra reported by Markevitc h et alJ 
(2002) for ObsIDs 3013 and 3419, which observed regions 
of the sky removed from any bright Galactic structures (see 
their Fig. 14). Those spectra were fit with a MEKAL ther- 
mal model, in both cases yielding a temperature of T ~ 0.2 
keV. We obtain the same result when fitting that model to our 
bright-September data. 

4.2. Predicted Intensities for Average Solar Wind 

To compare predicted absolute line intensities with ob- 
served values, we return to Eq. 1, and write the emissiv- 
ity as e,7 = v c n n n p fofiynO'i, where the ion density, is ex- 
pressed as the solar-wind (proton) density, n p , times the rela- 
tive abundance of oxygen, /<?, times the ion abundance rel- 
ative to O, /;. As a baseline for comparison with specific 
Chandra observations, we use average solar-wind parameters 
in the calculations that follow. As noted before, the average 
slow-wind velocity v c is 400 km s" 1 . Solar wind density is 7 
cm" 3 at 1 AU, the fractional O abundance is about 5.6 x 10~ 4 
ISchwadron & Cravens! I2000l) . and relative ion abundances 
are listed in Table[5] We assume that the ion density inside the 
magnetosphere is zero, and constant everywhere outside. Us- 
ing the undisturbed wind density is an approximation, as the 



wind sweeps around the Earth's magnetosphere in a wake- 
shaped structure (see Fig. 0, piling up on the leading edge 
of the bowshock, with higher-density lower-velocity shocked 
ions in the magnetosheath, but the resulting errors are compa- 
rable to other uncertainties in our model. 

The remaining factor in CT emissivity, and the one with 
the largest uncertainty, is th e neutral gas density. We use the 
analyt ical approximation of C ravens. Robertson. & Snowdenl 
(2001), based on Hodges' (1994) model of exospheric hy- 
drogen, which is that atomic H density falls off roughly as 
(l0R E /r)\ where R E is one Earth radius (6378 km) and r is 
the geocentric distance. On the leading edge of the magne- 
tosphere, at r ~ IQRe, n n is approximately 25 cm" 3 . On the 
flanks, the magnetosphere extends to r ~ \5Re- 

The brightness of a line I observed by Chandra while look- 
ing toward the Moon is given by 

1 [ DMmm -1 -2 -1 

Bji = — / eu(x)dx phot s cm sr (2) 
4?r Jo 

where x is the distance from Chandra. The neutral gas density, 
and therefore emissivity, is essentially zero at the distance of 
the Moon so we can replace Dmoou with oo in the integral. 
If we also approximate the look direction as being radially 
outward from Earth, then dx = dr and 
1 f°° 

Ba = 4^ J n n n pfofiV c yuVidr (3) 
_ vfipfofiyuajnM J°° ^ORe ^J 3 ^ ^ 

Vcn p fofiyn(7in„Q / 10R E \ 2 ... 
= 4^ 5Re \T~) ' ( ) 

"t/i \ I mm / 

where r m ,-„ is the geocentric distance to the edge of the mag- 
netosphere or the spacecraft's position, whichever is farther. 
The brightness of line / for average slow-wind parameters is 
then 

/ 1 A D \ 2 

B u = 9.95 x 10 14 yafm phot s _1 crrT 2 sr _1 . (6) 

\ 'min / 
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GSM X (R E ) 

FIG. 7. — Schematic of observation viewing geometry, with dotted curves 
for the July observations and solid for September. In Geocentric Solar Mag- 
netospheric (GSM) coordinates, Earth is at the origin, X points toward the 
Sun, and the XZ plane contains the Earth's magnetic axis. The Y and Z axes 
thus oscillate with a period of 24-hours, which makes Chandra's orbit appear 
non-elliptical. Positions of the magnetosph ere and bow shoc k are a pproxi- 
mate and based on the models of i Tsyganenko 1995, 1989) and I Bennett et al. 
1997), respectively. Chandra was barely inside the magnetopause in July, de- 
spite its apparent position in this projection. 



Chandra orbits between about 3Re and 2QRe, and the 
Moon observations were all made near apogee, so r mi „ = 
2QRe- {Chandra was in fact very slightly inside the magne- 
topause in July, but the difference in r, m „ is negligible.) Dur- 
ing both July and September, Chandra was looking toward 
the Moon through the forward flanks pointing slightly toward 
the leading edge (see Fig. [7}, and its orbit was inclined by 
~ 40°. A numerical integration taking into account the non- 
radial character of Chandra's, viewing angle results in a fac- 
tor of ~ 1 .5 increase over the radial approximation. A further 
slight adjust ment for the effects of the shock region (guided by 
the results of Robert son & Cravens! J200 3b)) brings the total 
increase to a factor of two, so that the final prediction for the 
Chandra observations, assuming a typical slow solar wind, is 
a line brightness of Bu = 5.0 x 10 yafiO'i phot s _1 cm" 2 sr _1 . 
The B mo dei values listed in Table0are the sums of Bu for lines 
within each line-group. 

It is often more convenient to compare predicted and ob- 
served counting rates rather than source brightnesses from 
spectral fits, particularly when energy resolution is limited, 
as with ROSAT. For an observation subtending a solid angle 
A57, the counting rate is Aft YluBuAu, where An is the in- 
strument effective area at the energy of line I and the sum is 
over all lines within the chosen energy range. For ACIS, one 
chip subtends 5.9 x 10~ 6 steradian, or 69.8 arcmin 2 , and we 
find that the predicted ACIS-S3 rate between 51 1 and 886 eV, 
encompassing all the He-like and H-like O lines, is 9.0 x 10~ 6 
cts s" 1 arcmin -2 . 

In comparison, the corresponding observed rate from the 
July S3 data, excluding ObsID 2469, is 34 ±20 x 10" 6 cts 
s" 1 arcmin -2 , while the ObsID 2469 rate is 155 ±49 x 10" 6 
cts s" 1 arcmin" 2 . The September rate, from the three bright- 
est ObsIDs, is 287 ± 39 x 10" 6 cts s" 1 arcmin" 2 , a factor of 



32 higher than predicted for an average solar wind, in gen- 
eral agreement, as o ne would expect, with the results from the 
spectral fits in 34.1H isted in Table[7] 

4.3. Adjustments for Actual Solar Wind 

An obvious question raised by the preceding analysis is 
whether or not the assumptions regarding solar wind parame- 
ters are appropriate for the July and September observations. 
This can be addressed by publically available data from solar 
monitoring instruments, specifically the Interplanetary Mon- 
itoring Platform 8 (IMP-8) 4 , and the Solar Wind Electon, 
Proton, and Alpha Monitor (SWEPAM) and Solar Wind Ion 
Composition Spectrometer (SWICS) onboard the Advanced 
Composition Explorer (ACE) 5 . Fig.[S]shows the available rele- 
vant data, namely the solar wind velocity and density, and var- 
ious element and ion relative abundances. Horizontal dotted 
lines denote the average solar-wind parameter values given 
previously. It is immediately obvious that CT emission rates 
should be much larger than average in September, and smaller 
in July. 

In July, the wind velocity v c averaged around 530 km s" 1 , 
versus 335 km s" 1 in September. Wind densities, on the other 
hand, are much higher in September than in July. The proton 
density measured by IMPS in July is ^3.5 cm" 3 , or half the 
baseline value, while in September the density ranges from 
~ 26 cm" 3 for the first ObsID (2468) to ~ 13 cm" 3 in the last 
(ObsID 3371). The AC£7SWEPAM densities are about half as 
large as the densities measured by IMPS, but are only avail- 
able in September, and only for level- 1 data which are not rec- 
ommended for scientific analyses. ACE also orbits around the 
Lagrange LI point lying 0.01 AU (~ 2AQR E ) toward the Sun, 
whereas IMPS is much closer to the Earth (orbital distance 
~ ?>5Re)- For all those reasons, we therefore rely on the IMP- 
8 densities. Note that ACE data plotted in Fig.[8]have been 
time-shifted to account for wind travel-time to Earth: July 
data by 1.584 x 10 6 km / 530 km s" 1 = 2990 s, and Septem- 
ber data by 1.435 x 10 6 km / 335 km s" 1 = 4280 s. IMPS 
data have at most a few hundred seconds of delay. After these 
corrections, the time-behavior of the IMPS and ACE density 
data match each other very well. 

Oxygen abundance data are not available for July, but the 
September data show a significant enhancement over average 
values, by roughly a factor of two during ObsID 2468 and ris- 
ing to more than a factor of three during ObsID 3371. In July, 
ACE/SWICS data indicate a relatively lower ionization state 
than usual for carbon (see Table This is also true of the 
September data, although the C 6+ /C 5+ abundance ratio is more 
volatile. There is a significant rise in the ratio during ObsID 
3371, which might be correlated with the increase in emission 
above 700 eV associated with especially highly-charged ions 
during that observation (see Table[3J. Overall, the abundance 
of the C 6+ ions that CT to produce C VI Lyman emission lines 
is lower than usual, in agreement with the relative weakness 
of those lines in the September spectra. 

AC£7SWICS does not measure the 8+ /0 7+ ratio, but the 
7+ /0 6+ ratio, which is usually 20:73 (see Table is less 
than half that in July and much higher in September. Deter- 
mining what the absolute ion fractions are can not be done 
precisely since the wind ions are not reall y in ionization 
equilibrium, but based on the tabulations of iMazzotta et alJ 
( 1998) we estimate that the 6+ :0 7+ :0 8+ relative abundances 

4 Data available from ftp://space.mit.edu/pub/plasma/imp/www/imp.html. 

5 Data available from http://www.srl.caltech.edu/ACE/. 
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are roughly 90:10:0 in July and 35:50:15 during most of the 
September observations, with a drop in ionization level to 
60:30:10 during most of ObsID 3370. That drop in the abun- 
dance of highly charged oxygen may be largely responsible 
for the lower X-ray emission observed during ObsID 3370 
(see Tables- 
Putting all the above factors together (v c , n p , rio/n p , and 
fi) and assuming that the July oxygen abundance was nor- 
mal, the CT X-ray intensity between 500 and 900 eV in July 
should be 0.25 times the baseline intensity, or 2.3 x 10~ 6 cts 
s" 1 arcmin" 2 . The predicted bright-September rate is 15 times 
the baseline rate and 60 times the July rate, or 140 x 10~ 6 cts 
s" 1 arcmin" 2 . Considering the model uncertainties, particu- 
larly our approximation for the neutral H density distribution, 
this is very good agreement with the observed quiet- July and 
bright-September rates of 34 ± 20 and 287 ± 39 x 10" 6 cts s" 1 
arcmin -2 , which are quoted with lcr uncertainties. 

4.4. Application to ROSAT Moon Observation 

We can also compare predicted and observed rates for the 
ROSAT Moon data. Because ROSAT was in low-Earth or- 
bit, observing perpendicular to the Sun-Earth axis through the 
flanks of the magnetosphere, r m ,„ = 15Re- (Note that the line 
brightnesses plotted in Fig. [5] are for this case, with the as- 
sumption of average solar wind parameters.) ACE had not yet 
been launched, but during the time of the observation, 1990 
June 29 2: 1 1-2:45 UT, IMPS measured v c = 400 km s" 1 with 
a slightly elevated proton density of 8.5 cm" 3 . Bu is then equal 
to 5.37 x 10 I4 y,v/;cr; phot s _1 cm" 2 sr _1 . 

Using the ROSAT PSPC effective area for A a , and with 
Ail = 3.0 x 10~ 5 steradians for the half-Moon, the observed 
rate, Ail Y^u^n^u, * s then predicted to be 0.0038 cts s , vs 
the vignet ting-corrected dar k-side observed value of ~ 0.20 
cts s" 1 ( Schmitt et a a difference of roughly a factor 

of 50. In the units of the ROSAT All-Sky Survey, the predicted 
CT rate is 1 1 x 10~ 6 cts s" 1 arcmin" 2 , while the observed dark- 
side rate is nearly 600 x 10~ 6 cts s" 1 arcmin" 2 . 

Given this l arge d iscrepancy, we look more closely at the 
ISchmitt et alJ {l991) observation. In their Fig. 4 we see that 
the intensity of full-band SXRB emission surrounding the 
Moon was ~ 2000 x 10~ 6 cts s" 1 arcmin" 2 , which is twice 
as high as the total rate of ~ 1000 x 10" 6 cts s" 1 arcmin" 2 
listed for the field of view in RASS map s during the ob serva- 
tion (RA, Dec - 1 l''44 m 30\ -03°22"). ISnowden et alJJl995h 
also noted this and concluded that "the lunar observation oc- 
curred during the time of a strong LTE." Although IMPS data 
on proton velocity and density do not indicate anything out 
of the ordinary at that time, we have seen from ACE data 
during the September 2001 Chandra observations that other 
wind parameters such as oxygen abundance and ionization 
level can have a very large effect on net CT emission. We 
therefore agree that the ROSAT data indicate an LTE, and 
deduce that the extra ~ 1000 x 10~ 6 cts s" 1 arcmin" 2 in the 
measured SXRB rate comprises ~ 600 x 10" 6 froma large in- 
crease in the geocoronal CT rate and ~ 400 x 10~ 6 from the 
associated transient excess of heliospheric CT emission. The 
rough equivalence of the geocoronal and heliospheric com- 
ponents of the LTE is consistent with model predictions by 
ICravens. Robertson. & SnowdenM2001l). 

5. IMPLICATIONS FOR THE SXRB 

5.1. Heliospheric Charge Transfer 

Although geocoronal CT accounts for much if not most of 
the X-ray emission during an LTE, and can be of roughly the 



same intensity as the cosmic X-ray background, its quiescent 
level is an order of magnitude or more below that of the typ- 
ical SXRB. Heliospheric CT emission, however, is several 
times stronger than qui e scent geocoronal emission, as was 
first shown by Cravens (2000), who presented a model for 
the distribution of neutral H and He within the heliosphere. 
Neutral gas is depleted near the Sun because of photoioniza- 
tion and charge transfer with H and He in the solar wind, and 
densities are higher upwind, with respect to the Sun's relative 
motion through the Local Interstellar Cloud, than downwind. 
Densities can be approximated by the relation n„ = «„oe~ A// ', 
where n„o is the asymptotic neutral density, r is distance from 
the Sun, and A is the depletion scale. 

Outside the heliosphere, in the undisturbed interstellar 
medium, the neutral H density is 0.20 cm" 3 , but behind the 
shock front the dens ity drops by nearly a factor of two, so 
that n„o = 0.12 cm" " 3 (C,loec k^re7aTin 997). Helium, in con- 
trast, is largely unaffected by the shock, and n„ o = 0.015 cm" 3 . 
Based on the theoretical work of Zank ( 1999), A for H is ap- 
proximately 5 AU upwind, 7 AU on the heliosphere flanks, 
and very roughly 20 AU downwind. He is less depleted 
near the Sun and has less sp atial variation, with A ~ 1 AU. 
Robertson & Cr avens! ( 1200 3a) have developed a more sophis- 
ticated model for the distribution of neutral H and He in all di- 
rections, but our predictions of absolute CT emission should 
not be significantly less accurate; in both cases, the estimated 
uncertainty is roughly a factor of two or three. Our calcula- 
tions, however, keep track of each CT emission line and the 
ROSAT PSPC effective area for that line energy, which will 
allow us to draw more specific conclusions regarding the rel- 
ative strength of emission in various energy bands. 

For heliospheric CT emission, Eq. 3 thus becomes 

1 f°° ( 1 AU\ 2 

Bu = —npfofiV c yuo-in n o e~ x/r dr, (7) 

4tt ViauV r J 

where we again approximate the look direction as radially 
outward from the Sun. Proton density, n p , decreases as the 
solar wind expands away from the Sun, leading to the r~ 2 fac- 
tor in the integral. If we isolate the ion-specific terms, yufiO~i 
(and ignore changes in the abundance of each ion species with 
distance from the Sun — ions change charge with each CT col- 
lision, but the path length for CT is many tens of AU), and 
evaluate the rest of Eq. 7 using the average solar wind and 
neutral gas parameters listed previously, we obtain 

Bit = Cyufidi phots" 1 cm~ 2 sr~ 1 , (8) 

where C has units of s" 1 cm" 4 sr" 1 and parametrizes the so- 
lar wind density and velocity and the neutral gas distribution 
along the line of sight. For CT with He, C = 1.76 x 10 15 , 
and for H it equals 4.37 x 10 15 , 3.12 x 10 15 , and 1.06 x 10 15 
when looking upwind, on the flanks, and downwind, respec- 
tively. Observations that look through the helium focusing 
cone downwind of the Sun, where He density is enhanced by 
roughly a factor of four (Michaels et al. 2002), will see more 
emission, with C ~ 10 16 . Note, however, that cross sections 
for He CT with the most important ions are roughly a factor 
of two smaller than for H CT (see Table |5| which also lists 
values of /; for the slow wind). 

Predicted ROSAT counting rates from quiescent CT emis- 
sion are listed in Table [8] The total of geocoronal and helio- 
spheric emission is around 80 x 10" 6 cts s" 1 arcmin" 2 , vary- 
ing by ±30 x 10~ 6 cts s" 1 arcmin" 2 depending on look direc- 
tion through the slow solar wind. The highest rate, roughly 
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FIG. 8. — Solar wind data for Chandra observations. Gray shading marks the duration of each ObsID. Dotted horizontal lines mark average values for the 
slow solar wind. Level 1 data from ACE (in September proton density plot) are not considered reliable, but are shown for comparison with IMPS data. IMPS 
and ACE proton velocity data in September are essentially coincident. Based on Eq. 1, and the corresponding data plotted here, the September X-ray rates are 
predicted to be much higher than in July, in agreement with observation. 
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TABLE 8. 
Model CT Rate Predictions 



Neutral 


Look 


Total ROSAT Rate 


Gas 


Direction 


(10 -6 cts s" 1 arcmin -2 ) 


Helio H 


upwind 


87 


Helio H 


flanks 


62 


Helio H 


downwind 


21 


Helio He 


any 


16 


Helio He 


He cone 


~90 


Geo H 


upwind 


24 


Geo H 


flanks 


11 


Geo H 


downwind 


~ 


Note. — 


- Results are for look directions within the slow 



solar wind. ROSAT always looked through the flanks of 
the Earth's magnetosphere; geocoronal rates for upwind and 
downwind directions are provided for comparison and appli- 
cation to other missions. Heliospheric emission when looking 
primarily through the fast wind is much lower. 



120 x 1(T 6 cts s" 1 arcmin , will be observed when look- 
ing through the He cone. For observations primarily through 
the fast wind, the geocoronal component remains the same, 
while the He and especially H heliospheric components will 
be significantly reduced, yielding a total of ~ 25 x 10" 6 cts 
s" 1 arcmin" 2 . For comparison, typical RASS rates away from 
bright Galactic structures are 600, 400, and 100 xlO" 6 cts 
s" 1 arcmin" 2 for the full, R12, and R45 energy bands, re- 
spectively. For lines of sight mostly through the slow wind 
(which is most of the sky during solar maximum, when the 
ROSAT background maps were conducted), our model there- 
fore predicts that combined quiescent geocoronal and helio- 
spheric CX emission accounts for roughly 13% of the total 
SXRB measured by ROSAT, 1 1% of the R12-band emission, 
and one-third of the R45-band emission. Overall uncertainty 
in our predictions is roughly a factor of two. 

Slightly more than half of the total CT emission (~ 55%) is 
predicted to be in the R12 band, which is consistent with the 
observatio n that the ROSAT LT Es appear most strongly in the 
R12 band iSnowden et al] 1 19951) . The true R12-band emis- 
sion is almost certainly even higher, however, because of the 
known incompleteness of our CT spectral model. Although 
we model Li-like O VI and Ne VIII emission, there is addi- 
tional L-shell (n — » 2) CT emission, mostly in the R12 band, 
resulting from Li-like and lower charge states of Mg, Si, S, Ar, 
and Fe. Although more difficult to model, these heavier ele- 
ments have significant abunda nce (in sum, ap proximat ely th e 
same as for bare and H-like C ( Schwadron & Cravens 2000)), 
and they have been invoked to explain much of the low-energy 
emission seen in comets ( Krasnopolskv 2004). 

Also note that some of the SXRB recorded in the RASS 
comes f rom point source s whic h could not be resolved by 
flOSAr- lMarkevitch et all fe002l) estimate 20%— so the frac- 
tion of the "true", i.e., diffuse, SXRB arising from CT is cor- 
respondingly higher. Whatever the average level of CT emis- 
sion is, we expect that it may well be the dominant contributor 
to the SXRB in some low-intensity fields, particularly in the 
R45 (3/4-keV) band. If so, CT emission will have a signifi- 
cant impact on our understanding of the SXRB and models of 
the local interstellar medium, particularly the Local Bubble. 



5.2. Once and Future Data 

ROSAT, of course, is not the only X-ray satellite affected 
by CT emission; any X-ray observation made from within 
the solar system will be impacted. Variations in the strength 
of background O emission have been noted (tentatively) in 
repeated Chandra observations of MBM 12, a nearby dark 
molecular cloud that shadows much of the SXRB (R. J. Edgar 
et al., in preparation), and in XMM-Newton observations of 
the Hubble Deep Field North, which also showed variable Ne 
and Mg emission (S. L. Snowden, 2004, private communica- 
tion). Such variations may be important when observing weak 
extended sources such as clusters where background determi- 
nations are critical. 

With future missions having much better energy resolu- 
tion, such as ASTRO-E2 with its 6-eV-F WHM-resolution mi- 
crocalorimeter array iStahle et alJ l2003). one should be able 
to see differences in background spectra between areas of 
the sky corresponding to low and high solar latitudes (around 
solar minimum) because of the difference in ionization lev- 
els between the slow and fast solar winds. It should also be 
easy to identify spectral features indicative of CT, such as an 
anomalously strong O VII Ka forbidden line (\s2s 3 S] — > Is 2 
^o). When excited by CT, that line is several times stronger 
than the resonance line {ls2p P] — > Is 2 l Sp), completely un - 
like their ratio in thermal plasmas ( Kharchenko et al. 2003). 
Other signatures of CT are the high-n Lyman transitions of 
O and C, which are strongly enhanced relative to Lya. In- 
deed, there are hints of high-n O L yman lines in a short (10 
sec) microcalorimeter observation (McCammon et al. 2002). 
There is also a strong feature at ~ 67.4 A (184 eV) in the 
Diffus e X-Ray Spect rometer (DXS) spectrum of the SXRB 
(Sand ers" et alJ lTOQl) which cannot be reproduced by ther- 
mal plasma models. A promising explanation is CT emission 
at 65.89 and 67.79 A (188.16 and 182.89 eV) from Li-like 
Ne VIII ls 2 5d -> ls 2 2p and H-like O VIII 5d, 5s -> 2p transi- 
tions, which are the strongest CT features within the 148-284- 
eV DXS energy range (see Fig. [5] and note that the plotted 
resolution is very similar to that for ASTRO-E2 and DXS). 

Although variability in foreground emission, whether from 
changes in the solar wind or in the position of the observ- 
ing spacecraft with respect to the magnetosphere and helio- 
sphere, can be an annoying complication when analyzing ex- 
tended sources, variations of CT emission in time and space 
also provide opportunities to learn more about the solar wind, 
geocorona, and heliosphere. Data from Chandra and XMM- 
Newton will provide important constraints on our models of 
CT emission, and future missions such as ASTRO-E2 should 
permit definitive tests. Given the unavoidable and rather large 
uncertainties in current theoretical models, observational data 
will be critical in this endeavor. 

6. CONCLUSIONS 

As described in this paper, we have detected significant 
time-variable soft X-ray emission in Chandra observations 
of the dark side of the Moon which is well explained by our 
model of geocoronal charge transfer. The observed bright- 
ness ranged from a maximum of ~ 2 x 10~ 6 phot s" 1 arcmin" 2 
cm" 2 , with most of the emission between 500 and 900 eV, to 
a minimum at least an order of magnitude lower. Predicted 
intensities, which are based in part on detailed solar wind 
data, match observation to within a factor of two, which is 
within the model uncertainty. Emission from O VII Ka, O VIII 
Lya, and a blend of high-n O VIII Lyman lines is detected 
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with high confidence, as well as probably Mg XI Ka and per- 
haps high-n emission from C VI. We also include estimates 
of heliospheric emission and find that the total charge trans- 
fer emission amounts to a substantial fraction of the soft X- 
ray background, roughly one-third of the rate measured in the 
ROSAT R45 band. Future observations with microcalorime- 
ter detectors should allow much more accurate assessments 
of the contribution of charge transfer emission to the SXRB 
because of its unique spectral signatures. 
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